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Abstract: The current status and prospects of deducing the proton-air cross section from fluorescence 
telescope measurements of extensive air showers are discussed. As it is not possible to observe the point 
of first interaction, Xi , directly, other observables closely linked to X\ must be inferred from the mea- 
sured longitudinal profiles. This introduces a dependence on the models used to describe the shower de- 
velopment. Systematic uncertainties arising from this model dependence, from the reconstruction method 
itself and from a possible non-proton contamination of the selected shower sample are discussed. 
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Introduction 

Indirect cosmic ray measurements by means of ex- 
tensive air shower (EAS) observations are difficult 
to interpret. Models needed for a deeper under- 
standing of the data have to be extrapolated over 
many decades in energy. This is the case for high 
energy (HE) interaction models, but also applies 
to the primary composition of cosmic rays. Un- 
fortunately a changing primary composition and 
changes in the HE interaction characteristics can 
have similar effects on EAS development and are 
difficult to separate. 

One of the key parameters for EAS development is 
the cross section <T p _ a ir of a primary proton in the 
atmosphere. Of course, only the part of the cross 
section leading to secondary particle production is 
relevant for EAS development, which we call for 
simplicity here (7 p - a ir- But also the production 
cross section contains contributions which cannot 
be observed in EAS. As diffractive interactions of 
primary particles with air nuclei do not (target dis- 
sociation) or weakly (projectile dissociation) influ- 
ence the resulting EAS, any measurement based on 
EAS is insensitive to these interactions. Therefore, 
we define an effective cross section to require an 



inelasticity fcj no i = 1 — 



of at least 0.05 



In the following the amount of traversed matter be- 
fore an interaction with fc; ne i > 0.05 is called X\. 
Taking this into account the reconstructed value 
of crp_ a i r needs to be altered by a model depen- 
dent correction a™^(k- me i < 0.05). This cor- 
rection amounts to 2.4 % for SIBYLL [1], 3.9 % 
for QGSJETII.3 [2] and 5.5 % for QGSJET01 [3], 
resulting in a model uncertainty of ^3 %. 

All EAS simulations are performed in the CONEX 
[4] framework. To account for the limited re- 
construction accuracy of a realistic EAS detector, 
I max is folded with an Gaussian function having 
20 gcm~ 2 width, which corresponds roughly to the 
resolution of the Pierre Auger Observatory [5]. 

X max -distribution ansatz 

The most prominent source of shower fluctuations 
is the interaction path length of the primary particle 
in the atmosphere. However the EAS development 
itself adds a comparable amount of fluctuations to 
observables like X max . This is mainly due to the 
shower startup phase, where the EAS cascade is 
dominated by just a few particles. Our approach 
to fit the full distribution of X max does therefore 
handle the primary interaction point explicitly and 



CpL air = 0- p _ air (fci ncl > 0.05). 
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Figure 1: Impact of a 10 % change of cr p _ a ir in 
QGSJETII at 10 EeV. Data from [6, 8, 9, 10, 11, 
12, 13]. 

the EAS development in a parametric way 



dP 



A* 



p— air 

Pax (AX + X s hift, A*„ Qir ) 



(2) 



where AX was introduced as X max — X%, Thus 
the X max -distribution is written as a double con- 
volution, with the first convolution taking care of 
the EAS development and the second convolution 
handling the detector resolution. In this model we 
have two free parameters A* _ air , which is directly 
related to cr*_ air , and X s hift, needed to reduce the 
model dependence. Note that Eq. (2) differs from 
the HiRes approach [6] and that used in the simula- 
tion studies in [7] by explicitly including the cross 
section dependence in Pax ■ 
The simulated Pax -distributions can be 
parametrized efficiently with the Moyal func- 
tion 



Pax (AX) = 



e 2 



(t+e-*) 



and t = 



AX 



/?V2tt 

using the two free parameters a and (3. 



ft 



(3) 
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Figure 2: Example fits of Eq. (3) to simulated 
Pax -distributions at 10 EeV. 

CONEX for several HE models such that the cross 
section used in the simulation is replaced by 

a— d (E) = a p _ air (P) • (1 + f(E)) , (4) 

with the energy dependent factor f(E), which is 
equal to for E < 1 PeV and 



f(E) = (/lOEeV - 1) 



log 10 (P/lPeV) 
log 10 (l EeV/1 PeV) 



(5) 



for E > 1 PeV, reaching f 1QEcV at E = 10 EeV. 
This modification accounts for the increasing un- 
certainty of Cp-air for large energies (see Fig. 1). 
Below 1 PeV (Tevatron energy), er p _ a i r is pre- 
dicted within a given HE model by fits to the mea- 
sured pp cross section. 

The cross section dependence of Pax and the cor- 
responding parametrizations are shown in Fig. 2. 































■* ' / 






A'' 




QGSJETII 

QGSJET 








SIBYLL 
□ o- pajr (QGSJET) 








a o- (QGSJETII) 








a (SIBYLL) 




1 , , , 1 


, , i , . 





800 820 
a [gem 2 ] 



Impact of <T p _ air on EAS development 

To include the cross section dependence of Pax 
in a cross section analysis at 10 EeV, we modified 



Figure 3: Resulting <7 p _ a ii -dependence of the 
parametrized PAx-distribution. The markers de- 
note the location of the original HE model cross 
sections. 
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Figure 4: Sensitivity and HE model dependence of 
the <7 p _air reconstruction for a pure proton compo- 
sition at 10 EeV. 

At large AX, the simulated distributions are not 
perfectly reproduced by the parametrizations. This 
effect worsens for large cross sections, as can be 
observed from the increasing x 2 /ndf (see Fig. 2). 
Also the deviation of the Moyal function from the 
Pax -distribution depends on the HE model. It is 
biggest for QGSJETII and smallest for SIBYLL. 
Unfortunately this disagreement produces a sys- 
tematic overestimation of ~30 mb for the recon- 
structed er p _ai r . This is visible in all the follow- 
ing results and will be addressed in future work by 
making the parametrization more flexible. 
The dependence of a and (3 on cr p _ a h can be in- 
terpolated with a polynomial of 2nd degree. Fig. 3 
gives an overview of this interpolation in the a-f3 
plane. Obviously the Pax predicted by different 
HE model are not only a consequence of the dif- 
ferent model cross sections. 

Results 

Pure proton composition 

In Fig. 4 we show the reconstructed Cptf a ; r f° r s i m ~ 
ulated showers with modified high energy model 
cross section , <J™°^ cd . The original HE cross 
section < 1 ° d i i r fied - a™*?} = can be recon- 

JJ — dil LJ — all 

structed with a statistical uncertainty of mb, 
whereas the uncertainty caused by the HE mod- 
els is about ±50 mb. At smaller cross sections 
the reconstruction results in a slight overestima- 
tion (< 50 mb). But for larger cross sections there 
occurs a significant underestimation of the input 
cross section. This is mainly due to the worse de- 



Figure 5: Systematic caused by photon primaries 
at 10 EeV. 
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Figure 6: Systematic caused by helium primaries 
at 10 EeV. 

scription of Pax by the used Moyal function for 
large values of <7 p - a ir (see last section). 
Photon primaries 

Primary photons generate deeply penetrating 
showers. Even a small fraction of photon show- 
ers has a noticeable effect on the tail of the X max - 
distribution [7]. Fig. 5 demonstrates how much a 
few percent of photons could influence the recon- 
structed Up-air. The current limit on the photon 
flux is 2 % at 10 EeV [14]. Note that there is a 
clear trend of an increasing x 2 /ndf with increas- 
ing photon fraction, meaning the photon signal is 
not compatible with the proton model. 

Helium primaries 

On the contrary, helium induced EAS are very sim- 
ilar to proton showers. Therefore their impact on 
Cp-air is significant and very difficult to suppress, 
see Fig. 6. Interestingly, even for large helium con- 
tributions there is no degradation of the quality of 
the pure proton model fit (x 2 /ndf is flat). Thus it 
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Figure 7: Composition impact on X max at 10 EeV. 

is not possible in a simple way to distinguish be- 
tween a 25% proton / 75% helium mixture or just 
a pure proton composition with a cross section in- 
creased by about 150 mb. 

Outlook: Mixed primary composition 

Fluctuations and the mean value of the X max - 
distribution are frequently utilized to infer the com- 
position of primary cosmic rays [15]. It is well un- 
derstood how nuclei of different mass A produce 
shower maxima at different depth X max (^4) and 
how shower-to-shower fluctuations decrease with 
A (semi-superposition model). 
The relative change of the X max -distribution from 
a pure proton to a pure mass A primary composi- 
tion can be evaluated using CONEX. To fit X max - 
distributions we use the formula [16] 



dP 



dX n 



N ■ 



(- 



V2(X„ 



ak+3-«) 



(6) 



with four parameters AT, X pca k, 7 and S. The nor- 
malization constant N was not fitted, but set to re- 
produce the known number of events. Fig. 7 shows 
how the X max -distributions for proton, helium and 
iron primaries are positioned relative to each other 
for several HE models. This relative alignment can 
be utilized during er p _ a i r -fits to reduce the com- 
position dependence. The total mixed composition 
X max -distribution is then the weighted sum of the 
individual primaries 



dP 

J YTnix 
UvV max 



(^max) 



dP 



dXr, 



' yAj , X max ) 



where the weights lo.i are additional free parame- 
ters to be fitted together with Xjhift and A*_ air . 
The shape of d ^ p (A) for A > 1 is always as- 
sumed to change relative to the proton distribution. 

First studies indicate that the correlation between 
the reconstructed composition and the correspond- 
ing (Tp-air does not allow a measurement of the 
cross section. The situation is expected to be more 
promising if the parameter X s hift is fixed, however, 
the model dependence of the analysis will then be 
larger than shown here. 
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